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The Jordan- Wigner transformation is known as a powerful tool in con- 
densed matter theory, especially in the theory of low-dimensional quan- 
tum spin systems. The aim of this chapter is to review the application 
of the Jordan- Wigner fermionization technique for calculating dynamic 
quantities of low-dimensional quantum spin models. After a brief intro- 
duction of the Jordan- Wigner transformation for one-dimensional spin 
one-half systems and some of its extensions for higher dimensions and 
higher spin values we focus on the dynamic properties of several low- 
dimensional quantum spin models. We start from the famous s = 1/2 
XX chain. As a first step we recall well-known results for dynamics 
of the 2-spin-component fluctuation operator and then turn to the dy- 
namics of the dimer and trimer fluctuation operators. The dynamics 
of the trimer fluctuations involves both the two-fermion (one particle 
and one hole) and the four-fermion (two particles and two holes) excita- 
tions. We discuss some properties of the two-fermion and four-fermion 
excitation continua. The four-fermion dynamic quantities are of inter- 
mediate complexity between simple two-fermion (like the zz dynamic 
structure factor) and enormously complex multi-fermion (like the xx or 
xy dynamic structure factors) dynamic quantities. Further we discuss 
the effects of dimerization, anisotropy of XY interaction, and additional 
Dzyaloshinskii-Moriya interaction on various dynamic quantities. Finally 
we consider the dynamic transverse spin structure factor 5zz(k,tj) for 
the s = 1/2 XX model on a spatially anisotropic square lattice which 
allows one to trace a one-to-two-dimensional crossover in dynamic quan- 
tities. 
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1. Introduction (Spin Models, Dynamic Probes etc.) 

The subject of quantum magnetism dates back to 1920s. E. IsinJH sug- 
gested a simplest model of a magnet as a collection of N spins which 
may acquire two values cr = ±1 and interact with nearest neighbors on 
a lattice as ^ JcriCj and with an external magnetic field as —hj^'^i- To 
explain the properties of the model we have to calculate the partition func- 
tion Z = Trexp(— which yields the Helmholtz free energy per site 
/ = hnijv^oo {—TlnZ/N) (in what follows we set fee = 1 to simplify the 

notations). In one dimension the problem was solved by E. Ising. Later 

9 

L. Onsager solved the square-lattice Ising model and we know the solu- 
tion in two dimensioni^l. There is no solution of the Ising model in three 
dimensions until now. 

Another version of interspin interaction was suggested by P. A. M. Dirac 
and W. Heisenberg. The Heisenberg exchange interaction reads ^ Jai-Oj = 
^ J {crfcTj + o'fcTj + o-fcTj) where the Pauli matrices a — {cr'-^ , , (t^) are 
defined as 

Denoting the halves of the Pauli matrices as — a" /2 (in what follows we 
set ?i = 1 to simplify the notations) we consider the following Hamiltonian 

We note that the Hamiltonian of the anisotropic XYZ Heisenberg model 
([2]) covers in the limiting cases some specific models like the Ising model 
(J^ = ^ 0), the isotropic XY (or XX or XXQ) model {J^ = jv , ,r = 
0), the anisotropic XY model (J"= ^ , = 0), the isotropic {XXX) 
Heisenberg model ( ~ — J^), and the Heisenberg-Ising (XXZ) model 
(J^ = jv = j, = A J). 

Again we would like to calculate the partition function Z of the spin- 
1/2 model Unfortunately this task is very complicated even in one 
dimension. Due to H. Bethe we know how to find the eigenstates of the 
spin-1/2 linear chain Heisenberg modeffl The famous Bethe ansatz for the 
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wave function has the form 



V") = •."-r)s„i •••s„J TT ■•• T), 

l<7ii<. ..<nr<N 



a(ni, ...,nr)= ^ exp i ^ k-p^rij + ^ X! ^^■'Pj ' 



2 cot — = cot — — cot — 
2 2 2 



Nk, = 2Tr\ + J20v, (3) 



where the sum in the definition of coefficients a{ni, . . . ,nr), V G Sr, runs 
over all r! permutations of the labels {1,2,..., r}, Vj is the image of j under 
the permutation V. For further details see, e.g., Ref. 

Let us briefly recall the quantities of interest in the statistical mechanical 
studies of the spin models. As we have mentioned already the thermody- 
namic quantities like the entropy, the specific heat, the magnetization etc. 
follow from the partition function Z = {~^>^/'^) — Trexp (— 

the sum runs over all states A of the system with energy E\. Usually 
we are also interested in the equal-time spin correlation functions, e.g. 
{si ■ Sj), ((. . .)) = Tr (exp {—(3H) (. . .)) /Z; their nonzero limiting values, 
(e.g., lim|j_j|_^oo (si • Sj)) may indicate the existence of long-range order in 
the system. 

Within a linear response regime we add to the Hamiltonian Hq a small 
perturbation Ho ^ Hq ~ b{t)B, where the external field b{t) couples to 
the dynamical variable B, and observe a response of a dynamical vari- 
able A, (Ait)) - {A)o = J^^ dt'xAB{t - t')h{t') with XAB{t - t') = ie{t - 
t'){[A{t), B{t')])Q (here 9{x) is the Heaviside step function). The Fourier- 
transform of the dynamic susceptibility XAB(i — t'), ^Xab('^) + '^'^Xab{'^), 
is the quantity which can be measured experimentally. We note that the 
real and imaginary parts of the dynamic susceptibility are connected via 
the dispersion (or Kramers-Kronig) relation. On the other hand, the imag- 
inary part of the dynamic susceptibility can be expressed with the help 
of the fluctuation-dissipation theorem through another dynamic quantity, 
the dynamic structure factor. Thus, Saa{^) — dtexp(iti)t) (A{t)A) — 
23xAA(w)/(l-exp(-/3w)). 

Usually, the operator A is constructed from the local operator of the 
considered system An as follows: Ak = {1/Vn) J2n=i exp(ifcn)A„. We can 
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also rewrite the dynamic structure factor in the following forms 




(4) 



A, A' 



Sometimes it is convenient to make the following change in the first line 
in Eq. dl: A„{t) A„{t) - (A), A„+i{0) A„+;(0) - (A). We also note 
that in the zero-temperature limit T — (or /3 — > oo) the second line 
in Eq. (gD becomes simpler, SAA{k,uj) = 27r I (GS|Afc|A) |^ (5 - wa), 
uj\ = E\ — Eqs ■ 

In what follows we discuss mainly the dynamic properties of spin-1/2 
XY chains; just for this class of spin models application of the Jordan- 
Wigner fermionization approach is most fruitful. We notice here that re- 
cently it has been found that CS2C0CI4 is a good realization of the spin-1/2 
XX chaiiP and calculations of the dynamic quantities for the correspond- 
ing spin models might be important for the interpretation of the data from 
dynamic experimenti^. As an example of earlier studies we may mention 
dynamic experiments on the spin-1/2 XX chain compound PrClJ^. 

The rest of this chapter is organized as follows. At first we briefly in- 
troduce the Jordan- Wigner transformation (Sec. [2|) and concisely discuss 
some of its generalizations (Sec.|3|). Then we consider in detail the dynamic 
structure factors for the spin-1/2 isotropic XY chain in a transverse field 
distinguishing the quantities which probe two-fermion, four-fermion and 
many-fermion excitations (Sec. |4|). Next we examine the dynamics for two 
slightly more complicated chains: the dimerized isotropic XY chain (Sec.[5|) 
and the XY chains with the Dzyaloshinskii-Moriya interaction (Sec. [6|). 
The results obtained for one-dimensional XY spin models do not involve 
any approximation. This is not true in the two-dimensional case for which 
the Jordan- Wigner approach provides only approximate expressions for dy- 
namic quantities. We illustrate the Jordan- Wigner fermionization approach 
in two dimensions examining some dynamic quantities for the square-lattice 
spin-1/2 isotropic XY model (Sec. [7|). We end up with a brief summary 



(Sec. El). 
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Jordan- Wigner fermionization 5 

2. The Jordan- Wigner Transformation 

To be specific, we consider tlie one-dimensional spin s ~ 1/2 XX Z Heisen- 
berg chain with the Hamiltonian 

N N 

H = Y,J {slsl^, + 44+1 + A«+i) -hY^sl- (5) 

n— 1 n— 1 

we imply either periodic or open boundary conditions in Eq. ([5]). Here the 
spin operators sf satisfy the commutation relations sf,Sj — i(5y e^^^sj, 
£0/37 is the totally antisymmetric Levi-Civita tensor with e^y^ = 1. In par- 
ticular, [sf,sj] = etc. Obviously, s" can be viewed as the halves 
of the Pauli matrices ([T]). After introducing the spin raising and lowering 
operators (or the ladder operators) — sf^ ± isjj (sj^ = (s+ -I- s~) /2, 
*rt — i^n " ^n) /^O the Hamiltoniau ^ becomes 

71 ^ 

+ A [snSn - ^ ) (4+lSrT+l " ^ ) ) " H ( ^ ^ ) ' 



We note that the spin raising and lowering operators satisfy commutation 
relations of Fermi type at the same site, i.e. 

{s,T,s+} = l, {s,T,s,T} = = (7) 

and of Bose type at different sites 

K > 4J = K> Sm] = [4> Sm] = 0> n 7^ m. (8) 

We may use the Jordan- Wigner transformatiorP to introduce Fermi 
operators according to the following formulas 

ci=s^, c„ = (-2sj)(-2s^)...(-24_i)s;, n = 2,...,iV, (9) 

c\ = s+, ct = i~2sl) (-2s^) . . . (-24_0 4, n = 2, . . . , TV. (10) 

(Sometimes one can find in Eqs. ([9]), pO)) instead of — 2s^ the identical 
expressions 1 — 2s+s^ = cxp (±i7rs+s~).) Really, the operators introduced 
always satisfy the Fermi commutation relations 

{Cra, cj„} — ^nmi {Cn, Cm} — {cjj, cj„} — 0. (H) 
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(To check this one has to note that (— 2s^) = 1 and that s^s^ = — s^s^.) 
The inverse transformation to the one given by Eqs. dH), ^ reads 




±mJ2c}cj]cl n = 2,...,N. (13) 



Moreover, the Hamiltonian ([6|) in terms of the Fermi operators (jU, ((TO 
has the following form 



H 



(we use c]c]+i = (~-^*|) ^j+i ~ "^tc). In the case of periodic 

boundary conditions implied for the spin Hamiltonian ^ the transformed 
Hamiltonian (|14p obeys either periodic or antiperiodic boundary conditions 
depending on the parity of the number of fermions. However, in what fol- 
lows the calculated quantities in the thermodynamic limit iV — s- oo will 
be i nsen sitive to the boundary conditions implied (for further details see 
Ref.^. 

From Eq. p4p it becomes clear that the spin-1/2 isotropic XY chain in 
a transverse (z) magnetic field with the Hamiltonian 

n n 

in the Jordan- Wigner picture is represented by the Hamiltonian 

and therefore is an exactly solvable modeilJJJ^. Moreover, the XY ex- 
change interaction may be anisotropic; then the intersite interaction has 
the form 

TX X X I jy y y 
^ (cicn+1 - C„4_^i) + I (44+1 - C„C„+i) (17) 
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Jordan- Wigner fermionization 7 

with J = (J^ + J^)/2, 7 = {J'^ — jy)/2. We can also consider an additional 
intersite interaction, the so-called Dzyaloshinskii-Moriya interaction, which 
does not spoil a simple ferniionic HamiltoniarP^ 

D - -> Y {cicn+, + c„cl^,) . (18) 

Moreover, within the Jordan- Wigner fermionization app roach we can ex- 
amine rigorously some types of multi-spin interactions^, for example, 



„x z y _ y z X 



^ (c,\c«+2 + c„ct +2) • (19) 



4 

Within the frames of the Jordan- Wigner approach we can also general- 
ize simple spin- 1/2 XY chains assuming regularly alternating Hamiltonian 
parameterJISl qj- gome types of random Hamiltonian parametersA^ and still 
face exactly solvable models. 

On the other hand, as can be easily seen from Eq. (|14p the Ising interac- 
tion between z spin components leads to interacting spinless fermions and as 
a result the advantages of fermionization are less evident. (Obviously, we can 
split the interaction term in the spirit of the Hartree-Fock approximatiorli^, 
however, the resulting theory will be only an approximate one. On the other 
hand, in the low-energy limit we can bosoniz e the ferm ionic Hamiltonian 
obtaining an exact low-energy effective theorjffSlinHO] -y^g cannot exam- 
ine rigorously within the Jordan- Wigner fermionization approach the case 
of the next-nearest-neighbor interaction since 

n { 1 - 24+iC„+i j Cn+2 - C„ ( 1 - 2cl^^Cn+l) C,^gD) 



J, J, _i_ ). r' 

^n^n+2 ^ ^n''n+2 '^n 

It is worthwhile to note here that recently the Jordan- Wigner fermioniza- 
tion approach has been applied to the spin- 1/2 isotropic XY model on a 
diamond chaiiJ^, however, the authors of that paper apparently missed 
some interaction terms in the fermionic Hamiltonian and their statement 
about rigorous results for such a model is wrong. Finally we note that an 
external magnetic field directed along x or y axes has an enormously com- 
plicated form in the Jordan- Wigner picture. 



3. Generalization of the Jordan- Wigner Transformation 

The Jordan- Wigner fermionization is a powerful tool for the study of quan- 
tum spin chains. Since the late 1980s there w ere several a ttempts to extend 
this approach to two (and three) dimensioni^^I^SIMIlS] .^g^j ^^^^ 
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values s > 1/^2H27I1H1 Por a review on the two-dimensional Jordan- Wigner 
fermionization approach see also Ref. 

Bearing in mind the Jordan- Wigner transformation in one dimension as 
a guideline we consider in the two-dimensional case the following relation 
between spin s = 1/2 and Fermi operators 

d- = exp (-iaj) s7 , 4 = ^^xp (iaj) st , 
s7 = exp (iaj) dj, st = cxp (-ia^) dl, 

Here d, d^ are the Fermi operators, the operators defined according to 



(1211) commute at different sites if the c-number matrix satisfies the 
relation 

exp(ii3^) =-cxp(ii?j^). (22) 

There are many choices of the matrix which realize the two- 
dimensional Jordan- Wigner transformation. Following Y. R. Wang^ we 
use the Cartesian coordinates i = (ix,iy) to construct a complex number 
T^; = ix + Uy = Irjl exp (i arg(rj)) and then choose 

% = arg (tj - r-) =3 In (tj - r-) = 3 In - + i{jy - iy)) ■ (23) 

Indeed, for such a choice Eq. (1^^ is satisfied, exp (^i^jT^ = 

exp (^i arg(T|- - Tj)j = exp (^i (^arg(Tj - r^-) ± t^) ) = - exp (i^ij) ■ Another 

1941 

choice of the matrix has the following forrri^ 

% - vr (0 {tx - jx) (1 - <5.,,, J + S,^,,J{iy - jy) (1 - 5,;„,,J) ; (24) 

here 9{x) is the Heaviside step function (see also Ref. ^^). 

After performing the Jordan- Wigner transformation (f2T|) for the two- 
dimensional spin-1/2 XX Z Heisenberg Hamiltonian one gets 

^ = E (i^ (4 '^^p {} ("J - ) + ^?°^p - "j) ) 4 
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with 

- ar^ = dr ■ A{r), 
A{r) = Va, = - ^ '.'^ 4^.- (26) 

(we have used Eq. (j23|l for a^? (|2T|) ). We need further approximations to 
proceed with statistical mechanics calculations for the Hamiltonian (I25|) . 
Within the mean-field description one assumes did? {dpdf) — (s5.) + 
1/2 1/2. We expect such an approximation to be valid in the case of 
zero magnetic field. For the mean-field description in the case of nonzero 
magnetic field and an analysis of the magnetization processes in the spin 
system see Ref. We also note that a more sophisticated (self-consistent 
site-dependent) mean-field treatment has been suggested as welPS. 

After adopting the mean-field approach we face the problem of particles 
in a magnetic field with the flux per elementary plaquette $0 = tt. We may 
change the gauge preserving the flux per elementary plaquette to make the 
Hamiltonian more convenient for further calculations. For example, for a 
square lattice we have 

« = E (4"! - ^^4) + (4'ir - 5) (4",- - 5 

m 

'^ij:-iy;ix,iy+i ~ Ji^ + 1 Ay +2 Ay — Ji^ + l,iy;i^ + l,iy+l — J- (27) 

In the one-dimensional case when either vertical or horizontal bonds vanish 
the Hamiltonian (P7|) transforms into Eq. (|14p (with h = 0). 

Recently A. Kitaev has suggested a new exactly solvable two- 
dimensional quantum spin modet^. This is a spin-1/2 model on a honey- 
comb lattice with interactions between different components of neighboring 
spins along differently directed bonds. An alternative representation of the 
honeycomb lattice is a brick- wall lattice (see Fig. [T|). The Hamiltonian of 
the model reads 

^= E ('^i-'*.- /^l+i.^ + -^^s.^^i.^, + J34;4/+i) ; (28) 

^'-|-/— even 

j and I denote the column and row indices of the lattice. We discuss in what 
follows a fermionic representation for the Kitaev modeP^. Let us perform 
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Fig. 1. A honeycomb lattice (up) with its equivalent brick- wall lattice (down). The 
bonds Ji run from south-west to north-east, the bonds J2 run from south-east to north- 
west, the bonds J3 run from south to north. 



the Jordan- Wigner transformation 

y y i k<i i<j 

(compare with Eqs. (HU, As a result we find that 

Ji 



(29) 



Ji s 



Jsslisli^^ J3 (a]jaj,i - 
Next we introduce the foUowing operators 

Cj,/ = a] , + flj,;, djj = i (^a]_j - flj,;^ , j + Z = odd; 

Cj,/ = i (^a] , - fljj^ , dj,i = a'^.j + aj,i, j + / = even. 
In terms of these operators the Hamiltonian reads as follows 



i+iaj,i+i - ^ ) -(30) 



H 



E 



j-t-/— even 



.Jl , . J2 I • "^3 n ^ 



(31) 



(32) 
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Since Dj i — idj^dj^^i are good quantum numbers the Hamiltonian p2p 
corresponds to a model of spinless fermions with local static Z2 gauge fields. 
Thus, Eq. ([5^ explains a hidden simple structure of the spin model 

The generalizations of the Jordan- Wigner transformation for arbitrary 
spin values were discussed by several author fP*^-^^-'^^', however, these 
mappings have not yet provided a substantial break-through for difficult 
strongly correlated problems. 



4. Spin-1/2 Isotropic XY Chain in a Transverse Field: 
Dynamic Quantities 

We start with the simplest spin-1/2 XY model, the transverse 
XX chain, with the Hamiltonian ([T5]). After performing the Jordan- 
Wigner transformation we arrive at a tight-binding model for spinless 
fermions (jl6p and after performing the Fourier transformation, Ck = 
J2n=i ^^P (i^") Cn (fc = 2Tm/N if the number of fermions is odd 
or k ~ 2TT{n + l/2)/N if the number of fermions is even, n = —N/2, —N/2 + 
1, . . . , N/2-1 if N is even or n = -(TV- 1)/2, -(A^- l)/2-M, . . . , (iV- 1) /2 
if A^ is odd), the Hamiltonian (fT6|) becomes diagonal 

H = Y,Ak(^clck~^^ , Ak^n + Jcosk. (33) 

As it has been mentioned above, for the analytical calculations discussed 
below we may consider only periodic boundary conditions for the fermionic 
Hamiltonian (i.e. k — 2Tin/N in Eq. (|33| ). 



4.1. Two-fermion excitations 

W e begin with the transverse dynamic structure factor Szz{k,ijj) 
^6 6 \ 64\6o\ rpj^g calculation of the zz time-dependent spin correlation 
function is straightforward. After exploiting the Jordan- Wigner trans- 
formation we have «(*)<+;) - «)«+;> (4 (t)c„(t)c|^_^;C„+;) - 

(4cn)(ci+,c„+z). Here ct(i) = (l/Viv) exp (ifcn) 4 W and 4(t) = 

c^exp(iAfet). Next we have to use the Wick-Bloch-de Dominicis theorem, 

(cliCfescIaCfe*) = (4iC/c2>(43C/c4> " (4-1 4-3 ) (^^2 Cfc4> + {cl^Ck^){ck^cl^), and 
to calculate the elementary contractions introducing the Fermi function 
nu = 1/ (1 + exp (/3Afc)), {cl^Ck^) = ^ifca'^fci, (4i42) = 0. As a result, the 
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final expression for the zz time-dependent spin correlation function reads 



<itK+i) - i-'^'-r = ]^ E (-i - 

fcl ,/C2 



• exp (i (Afcj - AfcJ rifci (1 - n^J , 

^^^) = ^E«)— ^Et-^^- (34) 



n— 1 k 



Plugging Eq. ([M)) into Eq. ([J]) we get the desired transverse dynamic struc- 
ture factor 

S,.{k,u;) = ^exp(-ifcO / diexpM)(«(t) - (s^)) - (s^))) 



dfeirifci (1 - ?ifci+fe) (5 + Afc, - Afci+fe) 
2|jsin|cos(| -|-fc*)l 



where — tt < fc* < tt are the solutions of the equation cu = 
-2Jsm{k/2)sm{k/2 + k*). 

The zz dynamic structure factor psp is governed exclusively by a two- 
fermion (one particle and one hole) excitation continuum. The propertie s of 
the two-fermion excitation continuum were discussed by G. Miiller et aP^l 
we present these results briefly below. The boundaries of the two-fermion 
continuum in the plane wave-vector k - frequency uj (we assume w > 0, 
— TT < /c < tt) are determined by the equations 

ui = -Aki + Ak2, k = -ki + k2{ mod (27r)), Afe = + Jcosfc, (36) 

where — tt < ki < tt. Moreover, in the ground state we have to require in 
addition n^^ > and 1 — n^^ > 0, i.e. A^^^^ < and A^^ > 0. 

We start with the zero-temperature case. In this case the two-fermion 
excitation continuum exists as long as |f2| < \J\. Let us introduce the 
parameter a — arccos (|ri|/| J|) and the following characteristic lines in the 
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13 



k-uj plane 



^i(fc) 
\J\ 

^2(fc) 
\J\ 



= 2 



= 2 



. k . 

sm — sm 
2 

. k . 

sm — sm 
2 



^3(fc) 
Ul 



sm- 



(37) 
(38) 
(39) 



The two-fermion dynamic quantities in the ground state may have non-zero 
values only within a restricted region of the k-u plane with the lower bound- 
ary uji{k) — wi(fc) and the upper boundary w«(fc) = uj2{k) if |fc| < tt — 2a or 
w«(A;) = 1-1^3 (fc) if TT — 2a < |fc|. Obviously, the two-fermion dynamic quan- 
tities may have only three soft modes ko = {0,±2a}. Moreover, there is a 
middle boundary of the two-fermion excitation continuum LOrn{k) = 0J2{k) 
if TT — 2a < |fc| along which the two-fermion dynamic quantities exhibit a 
jump increasing their values by 2. Finally, the two-fermion dynamic quan- 
tities show one-dimensional square-root van Hove divergencies along the 
curve LOs{k) — LO-^ik). In Fig. [2] we display the characteristic lines ([57|) . 




Fig. 2. The two-fermion excitation continuum which governs the ground-state two- 
fermion dynamic quantities. \ J\ = 1, \Q\ = 0.1 (a), |0| = 0.9 (b). We show the lower 
boundaries (bold lines), the middle boundaries (dashed lines), the upper boundaries (thin 
lines) and the lines of potential singularities (dotted lines). 



P9p which give the boundaries of the two-fermion continuum and potential 
soft modes and singularities. 

As temperature increases the lower boundary is smeared-out and finally 
disappears, the upper boundary becomes uj'i{k) along which van Hove sin- 
gularities occur. In the high-temperature limit the two-fermion dynamic 
structure factor becomes fi- independent. 
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In Fig. [3] we display the transverse dynamic structure factor Szz{k,i^) 




Fig. 3. Szz{k,u!) (gray-scale plots) for the chain l|15|l with J = —1, Q = (a), Q = 0.3 
(b), n = 0.6 (c) at T = and at T ^ oo (d). 



psp at zero temperature (panels a, b, c) and in the high-temperature limit 
(panel d). 

There are other dynamic quantities which probe the two-fermion exci- 
tation continuum. Let us consider the dimer operator 

Dn = «+i + Slsl^:^ -> i (4cn+l - C„cl+i) . (40) 

That operator is related to a perturbation to the Hamiltonian (|15p which 
mimics dimerization, cos(7rn)I?„. The dynamics of fluctuations of the 
dimer operator can be measured experimentally: the corresponding dimer 
dynamic structure factor is relevant to ph onon-assisted optical absorption 
processes in magnetic-chain compounds^. 

The calculation of the time-dependent dimer-dimer correlation function 
repeats all steps discussed above while deriving (|34p and ends up with 



{Dn{t)Dn+i) - (D)' 



-y 



COS exp (-1 (fci - fc2) «) 



kl ,k2 



■ exp (i (Afe^ - AfeJ t) rifci (1 - nfej , 



(D 



1 ^ 1 



COS k tanh ^ ^ . (41) 
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Inserting Eq. into Eq. (0]) we get the dimer dynamic structure factor 

SDD{k,uj) = ^exp(-ifcO / dtexp{iLjt){{Dn{t) - {D)) [D^+i - {D))) 
1=1 

dfci cos^ ( fci + ^ ) (1 - Uk^+k) 5{i^ + - Afci+fe) 



E 



cos 



s2 (I + k*) Uk* (1 - Uk+k*) 



(42) 



2 I Jsin I cos (I + fc*) I 
We can also calculate the zD and Dz dynamic structure factors 

S2D{k,uj) = exp (^i^^ J dfciCOS^fci + 

■riki (1 - rifei+fc) S{u + Afci - Ak^+k) , 

SD.{k,uj) = {SMk,^)T ■ (43) 

In Fig. [4] we display the dynamic structure factors SDD{k,uj) p2)) and 




Fig. 4. SDD{k,(^) (a, b) and SzDik,u)) (multiplied by exp(— ifc/2)) (c, d) (gray-scale 
plots) for the chain lO with J = -1. SDD{k,u}): O = 0.3, T = (a), T ^ oo (b). 
5^i3(fc,w): = 0.3, T = (c), n = 0.3, T = 0.1 (d). 



SzD{k, uj) (1^5]) at zero temperature (panels a, c), at low temperature (panel 
d) and in the high-temperature limit (panel b). 
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Comparing Eqs. ([55]) . (UH), we immediately recognize that all these 
dynamic quantities are governed exclusively by the two-fermion excitation 
conti nuum (for other two-fermion dynamic quantities see below and also 
Refs. 13311) and therefore all of them exhibit generic properties inherent 
in the two-fermion dynamic quantities. However, they also exhibit some 
specific properties originating from additional factors in the integrands in 
Eqs. (|35l) . ([42]) . ([43|l (e.g. singularities may be suppressed etc., compare 
Fig. and Fig. [3^). In general, the dynamic structure factor governed by 
the two-fermion excitation continuum can be written in the following form 

SAB{k,u!) ^ / dkidk2CAB {ki,k2) 

•n-fci (1 - nk^) S (u) + Afe^ - Ak^) Sk+ki~k-2,Q, 

CUkiM) = 1, 

CDD(fcl,fe) = COS^ 



2 ki + k2 



CzD{ki,k2) = i(cxp(-i/ci) + cxp (1/1:2)), 

Cm{ki,k2) = {CMki,k2)T . (44) 

All these quantities show generic properties (spectral boundaries, soft 
modes, singularity structure) and specific properties controlled by 

CAB{kl,k2). 



4.2. Four-fermion excitations 

We proceed by considering more complicate dynamic quantities. Namely, 
consider a trimer operator^ 

rp — „x X . y y 



That operator enters as a perturbation to the Hamiltonian (|15p which mim- 
ics trimerization, e^^^ cos(27rn/3)T„. The dynamics of fluctuations of the 
trimer operator, although it can be analyzed rigorously, is less evident from 
the experimental point of view. Its importance, however, is justified as a 
quantity of intermediate complexity between the zz and the xx and xy 
dynamic quantities. 

The calculation of the time-dependent trimer-trimer correlation function 
contains no new conceptual ideas but is somewhat tedious. The final result 
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for the time-dependent trimer correlation function reads 

{T,,{t)T,,+i) - {Tf = ^Y1 4t(^i. ^2) exp (-i (fci - k^) I) 

fcl,fc2 

• exp (i ( Afci - Afe2 ) t) rifci (1 - J 
+ X! C'tt(^1' ^2,fc3,^4)exp (-i (fci + fc2 - fca - fc4) /) 

• exp (i (Afei + Afe2 - Afc3 - A^J t) Uk^Uk^ {1 - Uk^) {1 - nkj , 

1 ^ 

= 1^ E (^«) = + - 2C0C2 (46) 



with 



d^^{ki,k2) = (1 - 2co)^ COS^ (fci + k2) 

+4ci (1 - 2co) (^cos^ (^/ci + y ) + cos^ + ^2 

+4cJ (cos^ /ci + cos^ k2 

I 9 \ ^ k-\ + ko 9 o fci — fco 

8 (-C2 + + 2C0C2) cos^ + 8cl cos^ 



2 

^ / 9 fci 9 fca 

+4ci (1 - 2co - 4c2) ( cos^ y + cos^ — 

+4c2 - 8ci - Sci + 4c2 + 16coCi - 8coC2 + 16ciC2, (47) 

C^1^(fcl,fc2,fc3,fc4) 
1C • 2 fcl - ^2 . 2 ^3 - ^4 2 ^1 + ^2 + fcs + ^4 . „ . 

= 16 sm sm cos > (48) 



and Cp = (1/^) 12k (P^) "-fc- Obviously, the calculation of such an aver- 
age as {cl.^cl.^Ck3Ckicl.^cl.^CkrCks) according to the Wick-Bloch-de Dominicis 
theorem is rather complicated. Substituting (|46)) into Eq. (|4]) we obtain the 
following result for the trimer dynamic structure factor 



STT{k,u;) = ^exp(-ifcO / dt exp{iujt) {{Tn{t) - (T)) {T„+i - (T))) 
1=1 

= sPj.{k,uj) + s'^^lik.o^Q) 
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with 



(2) I (2) 

J — 7T 

S^r^l{k,U!) ^ ^-^ J dkidk2dk3C^rp^{ki,k2,k3,ki + k2 - k^ + k) 

■Uk^rik^ (1 - riks) (1 - nfcj+fe2_fe3+fe) 5 + Afe^ + A^^ - Afeg - Kk^+k2-k3+&'^) 

For further details see Ref. — . 

As can be seen from Eqs. I^^ . ([501) . (EJ) the trimcr dynamic struc- 
ture factor is governed both by the two-fermion (one particle and one 
hole) excitation continuum discussed above, the term Srprp{k,uj), and by 
the four-fermion (two particles and two holes) excitation continuum, the 
term s!^^{k,uj). The four-fermion excitation continuum is determined by 
the conditions 

w = -Ak, - Ak2 + Afc3 -I- Afc^, k = -ki -k2 + fcs + kii mod (27r)),(52) 

where — tt < fci, A:2, ^3 < and Afe = + J cosk. Moreover, in the ground 
state we have to require in addition > 0, rifc^ > 0, 1 — nk^ > 0, 1 — Uk^ > 
0, i.e. Ak, < 0, Afc, < 0, Afc3 > 0, Afc, > 0. 

We start with the zero-temperature case. The lower boundary is given 
by one of the following curves 

^l'\k) ^ . \k\ . f \k\ 
' ^ ^ 2 sin sin ' " ' ' 



|J| 2 V 2 

;p^(fc) ^ k f \k\ 
' 4 cos - cos a -I- V 



,(3) 



|J| 4 V 4 

(fc) / / Ifcl 

^ ^ 2 sin a -I- y sin 2a + 



(k) ^ . ( \k\\ . \k\ 
^ - 2 sin a - -L-l sin 2a - -L-l 



,(4) 



|j| V 2 y V 2 



|j| =-4sin — sin(^a-— j (53) 

depending on the value of fi, |ri| < | J| and the value of fc, tt < fc < tt as 
is shown in the left panel in Fig. [5] The boundary between the region i 
(where Lof* (k) is the lower boundary) and the region j (where lof^ik) is 
the lower boundary) (see the left panel in Fig. [5]) is given by the formula 
|fc| = lij{a) where Iyi{ol) = 4 arctan ^ ^tan a — V tan^ a — 3^ /3^, |fc| < 
27r/3; /13(a) = tt - a, 7r/2 < |fc| < 27r/3; /14(a) = 2a; haia) = 27r - 
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Fig. 5. The lower boundary u)i(k) (left panel) and the upper boundary u)u(k) (right 
panel) of the four-fermion excitation continuum in the plane wave- vector k - transverse 
field O, I J| = 1. 



4a; ^34(0;) = + cosa - 1/2, 27r/3 < |fc| < tt; hz{a) = 4a (for further 
details see Ref. 1^^. The four-fermion contribution to dynamic quantities 
may exhibit soft modes |fco| — {0, 27r — 4a, 2a, 4a}. Next we pass to the 
upper boundary of the four-fermion excitation continuum which is given by 
one of the foUowing curves 



,(2) 



(fc) 



Ul 



1^1 

— 4 cos — cos 
4 



— 4 cos ■ 



(54) 



(55) 



depending on the value of il, \^\ < \ J\ and the value of /c, tt < fc < tt as 
is shown in the right panel in Fig. [S] The boundary between the regions 
1 and 2 is given by the curve |fc| = 4a. In Fig. [S] we compare the ground- 
state two-fermion and four-fermion excitation continua for two values of the 
transverse field il. The four-fermion excitation continuum always contains 
the two-fermion excitation continuum. The lower boundaries may coincide 
(e.g. in the zero-field case the lower boundary is | J| sin for both continua, 
panel a in Fig. [6^) whereas the upper boundaries are different. 

Next we turn to the van Hove singularities inherent in the four-fermion 
dynamic quantities. Evidently, the quantity 



■6 {u — \J\ cos ki 



S{k,uj) = / dfci 

J —TT 

\J\ COS k2 H- I J| cos ks 



dk2 / dk^S{ki,k2,k-i,k) 

J — TT 

-|J|cos(fc + /ci+fc2-fc3))(56) 



may exhibit van Hove singularities characteristic to the three-dimensional 
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Fig. 6. Lower boundaries and upper boundaries of the two-fermion and four-fermion 
excitation continua for \ J\ = 1 and = (a) and Q = 0.3 (b) at T = 0. The two-fermion 
continuum is shown shaded. 



density of states. The lines of potential singularities are as follows 

|J| -^"^'^2' 

\J\ 4 
ioi^\k) ^ k 



(57) 



The four-fermion dynamic quantities may exhibit cusp singularities (akin to 
density-of-states effects in thee dimensions) along these curves. We illustrate 
potential singularities in the frequency profiles of S{k,uj) at different 
k in Fig. [T] 

For nonzero temperatures the lower boundary is smeared out and fi- 
nally disappears. The upper boundary is given by Eq. (I54p . In the high- 
temperature limit the properties of the four-fermion excitation continuum 
become fi-independent. 

After discussing some generic properties of the four-fermion dynamic 
quantities (inherent in any four-fermion dynamic quantity) we illustrate 
some specific properties conditioned by the function Cy1^(fci, k2, fcs, k^) (|48p . 
In Fig. [8] we display the trimer dynamic structure factor (|49p . The con- 
tributions of the two-fermion excitation continuum and the four-fermion 
excitation continuum to this quantity can be easily distinguished. 

We may formally introduce the polymer operator 



p(0 



XX I y y 



(58) 
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Fig. 7. S{k,Lu) ((56]l vs at fc = 2tt/3 (a) and fc = vr (b) for S{ki, k2, ks, k) = 1 (bold 
curves) and S{ki, k2, ks, k) = nfe^n^^ (l-n^g) [l - Uk-^+k^-k^+k) , T = for Q = 
(solid curves), £7 = 0.3 (long-dashed curves), Q = 0.6 (short-dashed curves), Q = 0.9 
(dotted curves). Vertical lines denote the values of cj^"'' (k), j = 1,2,3 H57| l. 




Fig. 8. STT(k,u}) l l49t (gray-scale plots) of the chain JTSj with J = —1, S7 = (a), 
n = 0.3 (b), n = 0.6 (c) at T = and at T ^ oo (d). 



(evidently v'ii^ — Dn and Vii^' ~ Tn)- Now the dynamic polymer structure 
factor S-p-p{k, uj) will involve 2m-fermion excitations with m — 1, 2, . . . , L 
These quantities are of moderate complexity in comparison with Sxx{k,ijj) 
and Sxy{k, uj) which are enormously complex (see below). We also note that 
in the limit I oo 



,(2) 



(59) 
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The last term in Eq. ([55]) is nonzero only if ^ 0. 

In passing, we note that the niultimagnon continua of quantum spin 
chains have been discussed recently on general ground by T. Barnes ■ . Ob- 
viously, there is an essential difference in comparison with our case, since 
the Jordan- Wigner fermions obey the Fermi statistics and this point has im- 
portant consequences for the four-fermion excitation continuum considered 
in some detail above. 

4.3. Many-fermion excitations 

We pass to dynamic structure factors which are governed by many-fermion 
excitations. Let us recall that according to the Jordan- Wigner transforma- 
tion we have 



i (l - 2c\c,) . . . (l - 2ctic„-i) ( 



1 + - + - + 

< = -l^^^n (60) 

where we have introduced the operators ip^ = c}^ ± Cm- Obviously the 
operators (p^ are linear combinations of the operators Ck, in terms of 
which the Hamiltonian is diagonal (sec Eq. p3|) ). 

Consider now the xx time-dependent spin correlation function 

■Vt^l ■ ■ ■ ^^-l^J-l'P^y^J'fit+lVj+l ■ ■ ■ ^t+n-l'Pj+n-lVt+n)- (61) 

It contains a product of 2(2j + n — 1) ip^ operators (in contrast to 
4(s|(i)s|^„) = {ip^ {t)ipj {t)ip^_^_^ipj_^_^) which contains the product of only 
four (p^ operators). Therefore the calculation of xx and xy dynamic quan- 
tities (which are governed by many-fermion excitations) is essentially more 
complicated. 

Exact analytical results for xx and xy dynamic q uantit ies are rather 
scarce. At the high-temperature limit T — > oo we that 

2.2 



4(s,"(t)s^\„) = <5„,ocos(m)exp [-\j^t 



4(s^(04+„> = -'^n,o sin (nt) exp [-\j''t' ] . (62) 
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At zero temperature the xx and xy time-dependent correlation functions 
are extremely simple only when \ > \J\ ■ Consider for example the case 
ft > \ J\ when the ground state is completely polarized |GSs) = 11^^=1 I in) 
(in spin language) or completely empty CfejGSc) = (in fermionic lan- 
guage). Owing to the simplicity of the ground state s+ |GSs) = cJ„|GSc) = 
(i/VTv) X;fcexp(ifcm)4|GSc), s^|GSs) = 0. Therefore 

^{s'-W'^+n) = (GS.|s-(i)4+JGS,) 
— exp (ifcn — i (fi + J cos k) t) , 



4(s,^(04+n) ^ -i(GS,|,sT(t)4^JGS,) - ~4i{s^t)s^+,,) 



(63) 



and as a result 



Sxx{k, uj) ~ iSxy{k, uj) — —S (lu — fl — J cos k) . 



(64) 



Many results at T = refer to the as ympto tic behavior of the xx or xy 
time-dependent spin correlation function j^^ l ^^ i From the paper by A. R. Its 
et al we know the long-time asymptotic behavior at nonzero temperatures 



exp(/(n,0)), ^>1, 



(Sj (Os,+„> - I ^2(.^+,4) ^ ^ < ^ 



ZTT 



dp\n + Jt smp\ In 



tanh 



P {ft — Jcosp) 



27r 



tanh 



(65) 



On the other hand , we can obtain the xx and xy dynamic quantities 
numericalljH^'^*^ ^ 4 8|49|50 |51J Consider the slightly more complicated in- 
homogeneous spin- 1/2 anisotropic XY chain in a transverse field with the 
Hamiltonian 



N 



El jxx X X 

N N 



jxy „x V 



JV^ „V „x 



jvv „y „i/ 

•^3 -^J^j+l) 



1 / 1 \ 

1=1 i,i=i ^ 
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where 



To diagonalize a form bilinear in Fermi operators like (|66p we perform the 
linear canonical transformation 

N N 

Vk=Y^ {gknCn + hknci) , vl = (dLcl + KnCn) ■ (68) 

n—1 n—1 

The resulting Hamiltonian reads as follows 

N 

Ku in] flu — ^ 



{'y/c','?!-} = {?7/c',7?fe"} = {»7fc',7?I."} = (69) 

if the coefficients gi^^ i hkn satisfy the set of equations 

(gfchfc)M = Afc(gfchfc), 
gfe = {gki ■■■ 9kN ) , hfc = ( /ifci . . . /ifcAT ) , M = ^ ^ . (70) 

Further it may be convenient to introduce the linear combinations <i>fc„ = 
gkn + hkn and ^'fcn = Qkn — hkn which enter the relations 

N 

fj = 4 + - ^ i'^pjvl + %jVp) , 



N 

t 

p=l 



We calculate the time-dependent spin correlation functions using the Wick- 
Bloch-de Dominicis theorem. For example, 

= {'PnVn){'Pn+m'Pn+m) 
-{'Pn{t)^n+ni){^nit)fn+m) + if^ (t^n+m) i^n it)'Pn+,n) ■ (72) 
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The r.h.s. of Eq. ([7^ may be compactly written as the Pfafhan of the 4x4 
antisymmetric matrix 



Pf 



\-{'Pnit)'Pn+m) - {'fin it)'Pn+m) - i'fit+m'Pn+m) / 



4(.<(t)<+™) 

(73) 



Similarly (see Eq. (j6ip ). for the more complicated xx time-dependent spin 
correlation function we have 

4«(i)<+™> 



= Pf 



(74) 



\~{iptit)ip++J ^{ip^{t)ip++J -{ip+it)ip++J ... / 



i.e. (sj(i)*j+n) can be written as a Pfaffian of a 2(2j/ + rt — 1) x 2(2j + n— 1) 
antisymmetric matrix. The elementary contractions involved in (|73p . (j74p 
read 



((^+(t)^+> - E {"^PlKmmp) + %,%rnF{-Ap)) , 
p=l 
N 

p=l 
N 



p=l 

N 

(^7(i)^„) ^-Y.i'^PjKmFiAp) + *;,*p„.F(-Ap)) , 
p=l 

FiA.,) ^ j^^^^y (75) 
1 + exp (pAp) 

It is worthwhile to recall some properties of the Pfafhans which are 
used for calculating them. In the first numerical studies the authors used 
the relation 

(PfA)^ = detA (76) 
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and computed numerically the determinants which gave the values of Pfaf- 
fians according to ((76|) . On the other hand, the Pfaffian may be computed 
directljH^ noting that 



PfCU^AU) =detUPfA 



(77) 



and that 



Pf 



V 






Rl2 





. 


.0\ 


Rl2 








. 


. 











R34 . 


. 








—R34 


. 


. 











. 


•0/ 



— Rt2R: 



12-0,34 • ■ 



(78) 



We use the approach described to calculate the xx and xy d yna mic 
structure factors for the spin-1/2 transverse XX chain numerically^. To 
estimate the quality of the numerical procedure we compare our numer- 
ical findings with exact analytical results in the high-temperature limit 
and with exact asymptotics for finite temperatures in Fig. (S] Knowing the 










lo-ln 










' s 










Fig. 9. Panel a: Time-dependence of tlie autocorrelation function (sj (i)sj), j = 51 
at infinite temperature obtained numerically (symbols) and analytically (see Eq. (|62j) 
(solid curves). H = 2, 1 (downward and upward triangles), Q = 0.5 (open circles), Q = 
0.1 (squares), £7 = (filled circles). Panel b: Time-dependence of the real part of the 
autocorrelation function (sj(t)sj), j = 51 at £7 = for various temperatures obtained 
numerically (symbols) in comparison with asymptotics (|65^ . /3 = 5, 1, 0.1, 0.00001 (from 
top to bottom). The exact analytical result for /3 = is also shown (the lowest curve). 
Evidently only the slopes of the asymptotics should be compared with the numerical 
results. 



time-dependent correlation functions we obtain the corresponding dynamic 
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j+n, 



structure factors according to 

Sxx{k,u!) = exp (— ifcn) 25R ( / dt exp (i + ie) 

Sxyik,u;)^ J2 expHfcn)2i3^ dt exp (i (w + ie) (s^"(i)sj+„|^) 

with e +0. In practice we consider chains of = 400 sites, take j = 
41, 51, 61, n up to 50 or up to 100, and set e = ... 0.001 ... 0.1 (see 
Ref. El). The results of our calculations for Sxx{k,uj) and Sxy{k,uj) are 
illustrated in Figs. [lOl EI 




- ■ 3 G J ' 




Fig. 10. Sxxik,uj) (gray-scale plots) for the chain ||15|I with J = —1, for f2 = 0.0001 
(a), n = 0.3 (b), Q = 0.6 (c) at low temperature /3 = 20 and for fl = 0.6 in the 
high-temperature limit 13 = (d). 



Let us recall, the transverse dynamic structure factor Szz{k,uj) probes 
two-particle excitations, i.e. it is governed by the excitations which are 
composed of two Jordan- Wigner spinless fermions. The two-fermion exci- 
tation continuum has a sharp upper frequency cutoff at which Szz{k,uj) 
may diverge. At T = it has also a sharp lower frequency cutoff which 
touches = at fco (soft modes). Szz{k,uj) is almost structureless (apart 
from upper boundary singularities) and exists only for |r2| < |J|. In the 
high-temperature limit T ^ oo Szz{k, uj) becomes fi- independent. All these 
features are nicely seen in Fig. [3l 
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- ri 1- 




1.6 
] J 



(LE 

n4 

OJ, 




I.E 

14 



5;(A:,Lj) (gray-scale plots) for the chain II15II with J = —1 for Q = 0.1 (a, b), 
n = 0.3 (c, d) at low temperature /3 = 20. Positive parts are shown in panels a and c, 
negative parts are shown in panels b and d. 



In contrast, the dynamic structure factors Sxxik^uj) and Sxy(k,ijj) are 
many-particle quantities in terms of the Jordan- Wigner spinless fermions. 
The frequency range of these quantities is not a priori restricted, however, 
in the low-temperature limit T ^ Sxx{k, oj) and Sxy{k, w) are rather small 
(but nonzero) outside the two-fermion excitation continuum. These quan- 
tities show washed-out excitation branches roughly following the boundary 
of the two-fermion excitation continuum. [Although the results presented 
in Figs, [inilll] refer to the case J < (the ferromagnetic sign of the XX ex- 
change interaction) the results for J > (the antiferromagnetic sign of the 
XX exchange interaction) follow by symmetry. In fact, while changing the 
sign of XX exchange interaction, -I- J — J, we get Sxxik^uj), Sxy{k,Lj) 
given by Eq. (|79p in which the wave- vector is changed fc ^ fc=p7r.] From the 
exact calculation in the strong-field zero-temperature limit ([M]) we know 
that Sxxik^Lo), Sxy{k,ijj) are proportional to 5 (w — A^), A^ = + Jcosfc. 
In the high-temperature limit T oo Sxx{k,uj) and Sxy{k,u}) become k- 
independent, but depend on f2. All the features described can be seen in 

Figs. [la [m 

It is instructive to compare our precise numerical findings in the low- 
temperature limit with the resul ts for the ground-state dynamic structure 
factors obtained by bosonizatiorffSH^I^. In the case £7 = within the 
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framework of the bosonization approach we have 

9{uj-\vk\) 



a 



X, z. 



(80) 



Here v = J is the velocity and rj^ = 1/2, rjz = "2, are the exponents which 
describe correctly the singularity at the lower continuum boundary wi (fc) = 
I J sin fc I \Jk\ as fc ^ or A: ^ ±7r. In the case of nonzero transverse fields 
f2 ^ the values of the Fermi momentum and Fermi velocity are changed. 
In Fig. [T^l we compare the predictions of the bosonization approach ([50]) 




Fig. 12. Sxx{k,u)) for the chain 1151 1 with J = —1; frequency profiles at A; = 
0, 0.1, 0.2, 0.3 (from left to right) at f2 = (panel a) and C = 0.3 (panel b). Bosoniza- 
tion results which follow from Eq. II80I I are shown by thin lines (v = 1 for = and 
V = 0.9539 . . . for S7 = 0.3); numerical results at low temperature j3 = 100 are shown by 
solid lines. 



with the numerical results at low temperatures. 

To summarize, in this section we have discussed dynamic properties of 
the spin-1/2 transverse XX chain within the Jordan- Wigner fermioniza- 
tion approach. Within this scheme the spin Hamiltonian corresponds to 
the Hamiltonian of noninteracting spinless fermions. The transverse dy- 
namic structure factor corresponds to the fermionic density dynamic struc- 
ture factor and probes the two-fermion excitation continuum. There are 
more dynamic structure factors which probe the two-fermion excitation 
continuum, e.g., the dimer dynamic structure factor. All two-fermion dy- 
namic quantities have common features (spectral boundaries, potential soft 
modes and singularities) and specific features. There are also dynamic quan- 
tities which probe the four-fermion excitation continuum; as an example of 
such a quantity we have discussed the trimer dynamic structure factor. 



December 27, 2008 16:39 



WSPC/Trim Size: 9in x 6in for Review Volume 



oderzh'070612 



30 O. Derzhko 

Remarkably, the dynamic structure factors which are associated with the 
dynamics of fluctuations of the x or y spin components (in contrast to 
the transverse dynamic structure factor which is associated with dynamics 
of fluctuations of the z spin component) are enormously complex within 
the Jordan- Wigner description since they probe many-fermion excitations. 
Nevertheless the two-fermion excitation continuum is still important for 
these dynamic quantities at low temperatures. As we have observed in our 
numerics, most of the spectral weight is concentrated along the boundaries 
of the two-fermion excitation continuum (it was also noted earlier for the 
XXZ Heisenberg chain'''^). This is not the case in the high-temperature 
limit when these dynamic structure factors show Gaussian ridges. In the 
next two sections we shall follow to what extent our observations survive 
for more complicated spin-1/2 XY chains. 

5. Dimerized Spin-1/2 Isotropic XY Chain in a Transverse 
Field 

Now we pass to the dimerized spin-1/2 XX chain in a transverse fleld. The 
Hamiltonian of the model reads 

= ^ J (1 - {-ITS) + + E < 

n n 

E ^ (1 - (-1)"'^) (c,\c„+i - c„cl+i) + f^E (cic^ - ^) , (81) 

n n ^ ^ 

where 5 is the dimerization parameter (0 < (5 < 1). After performing 
consequently the Fourier transformation, = (^1/\/N^ '^j^exp (ikn) c^, 
k = 2np/N, p = -N/2, . . . ,N/2 - 1 (TV is even), and the Bogolyubov 
transformation, Ck — Uk+T,iik + Wk^k+^m Uk — (1/V2) ^1 + |cosA:|/efc, 
Vk = sgn (sin(2A;)) (l/\/2) y^l — \ cosfc|/efe, ek = \J cos^ k ^ sin^ fc, the 



Hamiltonian becomes diagonal, H = '^k ^k {vl^k — with the ele- 

mentary excitationenergy Ak = ^ + Xk, A/c = sgn(cosfc)Jefe (for further 
details see Refs. I^^IMI). 

The calculation of the transverse dynamic structure factor follows the 
scheme described in Sec. 2] and ends up with the following result 

Szz{k,Lu) 

dfci [iuk^Uki+k + Vk^Vk^+kf nki (1 - rik^+k) 5 + Xki - ^ki+k) 

+ {ukiVk^+k - Vk^Uk^+kf nki (1 - Uk^+k+Tr) S{UJ + Afc, - Xk^^+k+TrJ^^) 
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(see alsoRefs.l3^ESE^. 

Again for the xx and xy dynamic structure factors exact analytical 
results are rather scarc e. In the high-temperature limit only the autocorre- 
lation functions survivj^ 



61 { J+i, ^ ) H, (j+t, -Ij- 



•exp I -im - I 1 - I J^t' I . (83) 

Here the Jacobian theta and eta functions are defined as follows 

CXD 

01 (u, k) = ^ exp (2niz) , 

n— — OQ 

oc 

Hi{u,k)= q("+^)%xp((2n + l)iz), 



n— — 00 



= 1^-" ;c(fc) ) ' ^=2^fc) 

and the complete elliptic integrals of the 1st and the 2nd kinds are given 
by 

/C(fc) = E{k) ^ r d9^/\-k-^ain^e. (85) 

^0 Vl - /c2 sin^ e Jo 

Moreover, = (1 ± ,5)^ /4. 

In the strong- field limit > | J| at T = we can repeat the calculation 
of the previous section to find, for example, for the xx time-dependent 
correlation function and the corresponding dynamic structure factor the 
following result 

4(s^(0s^+„) ^ Xl'^^P (^^'^) ("feexp (-iAfei) + vl exp (-iAfe+^t) 

k 

-i (-1)-'+" UkVk (exp (-iAfci) - exp (-iAfc+^i))^ , (86) 
S^^k, ^) = I {^l^ - ^k) + vis {lo - Afc+^)) . (87) 
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For arbitrary values of temperature and transverse field the xx and xy 
dynamic structure can be obtained numerically^^. 

Let us discuss the dynamic quantities of the dimerized transverse XX 
chain. We begin with the transverse dynamic structure factor which can be 
written as (compare with Eq. (|44|) ) 



S,,{k,uj)= / dA:idfc2CW(A:i,fc2) 

J — TT 

■riki (1 - rifej (5 (w + Afej - AfcJ 4+fei-fe2,o 
dkidk2C^^Hki,k2) 

■nki (1 - Uk^) 5 + ^ki - AfcJ 4+/ci-fe2+7r,0, 
C^^^ (fcl , fc2) = (Ufci + f 

C^'^\ki,k2) = {uk^Vk^ -VkiUk^f . 



As can be seen from Eq. (I88|) the transverse dynamic structure probes two- 
fermion excitations. Szzik,i-u) may have nonzero value within a restricted 
region of the k-uj plane when there is such a wave- vector ki, —n < ki < tt 
that uj — —Xki + Xki+k or uj = —Xki + Xki+k+Tr- Moreover, at zero tem- 
perature there are additional restrictions arising from the Fermi functions. 
The lines of potential singularities follow from the analysis of the equa- 
tions dAfej/dfci — dXki+k/dki = and dXk^/dki — dAfcj+/c-|_7r/dfci — 0. The 
characteristic lines in the k-uj plane which determine the behavior of the 
transverse dynamic struct ure factor were reported for the first time by 
J. H. Taylor and G. MiiUei'^. 

In Fig. [13] we show the region of the k-to plane in which the two- 



^^^^ '^j— lyr- 



II 







Fig. 13. Location of the roots of equations uj = — A^^ +^ki+k E'nd = ~^ki -l--^fci+fc+7r 
(— TT < fcl < tt) in the k—uj plane for 5 = (panel a) and 5 = 0.1 (panel b): light region: 
no roots, light gray region: two roots, dark gray region: four roots. 
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fermion dynamic quantity may have nonzero values. In Fig. [T3]we show the 




-.■i O -I n L 2 3 -.T -2 -I D I 2 3 




-3 -2 -I n I I 3 -3-2-10123 



Fig. 14. Szz(k,u)) (gray-scale plots) for the chain 18111 with J = —1, 5 = 0.1 at different 
temperatures /3 = oo (a, b), /3 = 20 (c, d), /3 = 1 (e, f) for f2 = (left panels a, c, e) 
and O = 0.3 (right panels b, d, f). Note that the results at /3 = 1 for Q = and SI = 0.3 
(panels e and f) are practically indistinguishable. 

transverse dynamic structure factor at different temperatures. Comparing 
Fig. 113b and Fig. [T3] one can see the van Hove singularities and the effects 
of the Fermi functions and the C*^^-*- and C^^-' -functions. In Figs. [TH [16] we 
show Szz{k^ Lo) at various values of the transverse field for two temperatures 
/3 = oo and /3 = 20. 

Next we pass to the xx dynamic structure factor obtained numerically. 
Typically we consider chains of = 400 sites assume in (jll) n = 41, / up 
to 50, consider t up to = 200 and take e = COOllSSl. In Fig. [17] we show 
Sxx{k,u}) of the dimerized transverse XX chain at low temperatures for 
different values of the transverse field. 
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In contrast to the zz dynamic structure factor which is a two-fermion 
dynamic quantity, the xx dynamic structure factor is a many-particle dy- 
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a 1 



□ 




-3-2-10113 -3-2-10123 




Fig. 17. Sxx{k,uj) (gray-scale plots) for the chain H81| l with J = —1, 5 = 0.1 at low 
temperature /3 = 20 and for Q = 0.1 (a), Q = 0.3 (b), f2 = 0.6 (c) and in the high- 
temperature limit (3 = for Q = 0.6 (d). 

namic quantity within the Jordan- Wigner fermionization approach. There- 
fore, nonzero values of Sxx{k,u)) far above the two-fermion excitation con- 
tinua may be expected. However, as can be seen in Fig. [17] the opposite is 
true: At low-temperatures Sxx{k,(-u) shows several well-defined excitation 
branches which follow roughly the boundaries of the two-fermion excita- 
tion continua. Although we can describe the low-energy physics also using 
the bosonization treatment, high-frequency features cannot be reproduced 
within such an approach. 

Finally we note that the dimerized XX chain does not show bound- 
state branches; within the fermionization approach this may be related to 
the absence of interactions between fermions. In contrast, a particle-hole 
bound state can be observed in the dimerized Heisenberg chairPZl. 



6. Spin-1/2 XY Chains with the Dzyaloshinskii-Moriya 
Interaction 

In this section we examine the effect of the Dzyaloshinskii-Moriya interac- 
tion (actually, the z component of the vector of the Dzyaloshinskii-Moriya 
interaction, see Eq. (fTSjl ). The Hamiltonian of the transverse XX chain 
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with the Dzyaloshinskii-Moriya interaction reads 

H = J2iJ + 'n^n+l) + D -hY.<- (89) 

n n 

Interestingly, the Dzyaloshinskii-Moriya interaction can be eliminated 
from the Hamiltonian (I89p resulting in renormalization of the isotropic XY 



exchange interaction''^. To see this, consider the following spin axes rotation 



Sn = Sn COS (pn 



s„sm0„ + si; cos0„, 



5„ ^ 5„ S„ 



0„ = (n - 1)(^, tan(^=y. (90) 
After such a unitary transformation the Hamiltonian (I89p becomes 



J = sgn( J) yJ,P + 1)2. (91) 



Using the unitary transformation (|90p we can examine the effect of 
the Dzyaloshinskii-Moriya interaction using the dynamic quantities of the 
transverse XX chain without the Dzyaloshinskii-Moriya interaction dis- 
cussed already in Sec. H) First of all we note that the zz dynamic structure 
factor is given by Eq. ([35|l . however, with = —h + Jcosfc. The for- 
mulas determining the two-fermion excitation continua boundaries are still 
given by Eqs. ([37]), (|38ll . p9|) but with J instead of J on the l.h.s. of these 
equations and in the definition of the parameter a. 

Exploiting (|90|) we find the following relations between the xx and xy 
dynamic structure factors of the model ((89|l (l.h.s. of Eq. ([921) ) and the 
dynamic structure factors of the model ([OT|) (r.h.s. of Eq. ([M)) )^ 

Sxx{k,uj) = ^ {Sxx{k - + Sxx{k + 

+ i {Sxy{k - ip,Uj)\j - Sxy{k + ip,U!)\j)) , 
Sxy{k,U!) = i {Sxy{k - + Sxy{k + <i5,Cj)|j 

-i[Sxx{k - - Sxx{k + (p,uj)\j)) ■ (92) 
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Therefore, using Eq. ([S^ we obtain for the model ([55]) 



4J 



exp 



{lu + hY 



V^/ / (c. + hY 

\Sxy[k,uj) = ^ I exp I — 



exp 



exp 



{lo - hf 

{uj - hY 



(93) 



i.e. in the high-temperature hmit Sxx{k, uj) and Sxy{k, oj) are fc-independent 
and display a single Gaussian ridge at oj = \h\. 

In the zero-temperature and strong-field limit (T = 0, \h\ > \/lP + D^) 
according to ([5^ and ([M)) we find 



Sxx{k,Lu) = -sgn{h)iSxyik,uj) 



= —6 [uj — \h\ — J cos {k + sgn{h)tp) 



(94) 



For arbitrary values of temperature and transverse field we use Eq. 
(|92p and numerical results for the xx and xy dynamic structure fac- 
tors of the transverse XX chain (PT|) (see Sec. 0]) to reveal the effect of 
the Dzyaloshinskii-Moriya interaction. Some of our findings are plotted in 
Fig. [T8l where we show Sxx{k,uj) for D = (left panels) and D ^ (right 




1 1 1 [ 1 r 



> ■ 4 




Fig. 18. Sxx{k, u>) (gray-scale plots) for the model II89I I with J = —1, D = (left panels 
a, c) and D = I (right panels b, d) for h = (upper panels a, b) and h = —0.6 (lower 
panels c, d) at low temperature (3 = 20. 
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panels) at different values of the transverse field h. 

We recall that in the low-temperature limit when J < and D — 
Sxxik, uj) and Sxy{k, to) are concentrated in the k-uj plane along the curves 
(f37|l . ([38|) . ([39|l which determine the boundaries of the two-fermion excita- 
tion continuum uji{k), ujm{k) and uJu{k) (see Sec. IH Figs. [101 E]). [For the 
antiferromagnetic sign of XX exchange interaction J > these dynamic 
quantities are concentrated along the curves uji{k ± tt), LOm{k ± tt) and 
a;„(fc±7r) as it follows from simple symmetry arguments.] In the case when 
the Dzyaloshinskii-Moriya interaction is present, D 0, the two-fermion 
excitation continuum splits into two continua (see Fig. [T5)) . the 'left' one 
with the boundaries uJi{k — ip), uJm{k — (p) and uju{k — ip) and the 'right' 
one with the boundaries uji{k + ip), LOm{k + ip) and w„(fc -|- ip). (The 'left' 
and the 'right' continua are connected by symmetry operation.) The larger 
D is the larger is the splitting controlled by <y9 = arctan(Z)/ J). At fixed 
D ^ and h — the spectral weight is equally distributed between the 
left and the right continua (panel b in Fig. [T8)) . While \h\ increases from 
to \/\P + D^ the spectral weight moves from one continuum to another 
continuum (panel d in Fig. llSp . 

We note in passing that the discussed peculiarities of the xx dynamic 
structure factor may be used for an unambiguous determination of the 
Dzyaloshinskii-Mor iya interacti on in chain compounds for example, in res- 
onance experimentJSnEnESISI] 

In the case of the anisotropic XY chain the Dzyaloshinskii-Moriya in- 
teraction cannot be eliminated by the transformation (P(7|) . Now we face the 
Hamiltonian 

n n 

E/ J ^iD ^ J — \D ^ 7/^1 \ 

n ^ 

+n (^ctc„-ijj(95) 

with J = (J^ -f jy)/2 and 7 = (J^ - jy)/2. This Hamiltonian can 
be put into a diagonal form by performing the Fourier transforma- 
tion, Cfe = (^l/VlV^ X]„6xp (ifcn) c„, c„ = I]fcexp(-ifcn)cfc, 
and the Bogolyubov transformation, Ck = —iukPk + Vk0^_k^ Pk = 
iufeCfe -f WfecLfe, Uk = sgn (7sinfc) (l/\/2) ^1 + (il + J cosfc) /Afc, Vk = 

(l/\/2) Vl - {^ + Jcosk) /Afc, Afc = J{n + Jcosk)^ +-1"^ sin^ fc. The fi- 
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nal result reads H = I]^, (^/S^/Sfc - 1/2^, Dsinfc + A^. We no- 

tice that the elementary excitation energy spectrum is gapless when il^ < 

J2 _|_ £)2 _ ^2 g^j^j ^2 < jj2 ^jjgjj ^2 ^ j2 g^j^^ ^2 ^ £)2 

The calculation of the transverse dynamic structure factor follows the 
lines explained in some detail in Sec. |4]and ends up witll^ 



3 

S'£)(fc,t^) = J dkiB''^\ki,k)C'-^\ki,k)S (^Lu - E''^\ki,k)y 
Bi^Hk,,k)^Bi^Hh,k) = WI^, Bi^Hh,k) = '-±1^1^, 

{n + J cos (fci - I)) (f2 + J cos (fci + I)) - T^sin (fci - |) sin (fci + |) 

C(^)(fcl,fc) = (^1 + 

C^^\ki,k) = n^^^_kn_f^^_k, 
i?«(fci,fc)-A,^+. +A_,^+., 
i?(2)(fci,/c) = A,^+. -A,^_., 
E^'\k^,k) = -A,^_.-A_,^_(0B) 

The transverse dynamic factor, as it follows from Eq. (j96p . is shown in panel 
c in Fig. [19] for a typical set of parameters. 

From Eq. (j96p we see that the transverse dynamic structure factor is 
governed by three two-fermion excitation continua. Let us discuss some 
properties of these continua (see Fig. [20] where we show two-fermion exci- 
tation continua for a specific set of parameters J = —1, 7 = 0.5, D = 1, 
fl = 0.5). We begin with the high-temperature limit when the Fermi fac- 
tors are not essential (left panels in Fig. I20p . The two-fermion dynamic 
structure factor may have nonzero values in the k-uj plane if the equation 
uj — E^J\ki,k) = has at least one solution ki, — tt < fc^ < tt. Next, the 
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Fig. 19. Sxx{k,Lij) (a), Syy{k,Lo) (b), Szz(k,ui) (c) for the spin chain II95I I with J = — 1, 
7 = 0.5, D = 1, f2 = 0.5 at low temperature /3 = 50. Note that these quantities are 
shown for k that varies from — tt to Stt. 

lower and the upper boundaries are given by 

Oj[^\k)^ min {Q,E^^\ki,k)\ , 

— 7r<ki<7T L J 

ujiPik)^ max {E^^^kuk)] . (97) 

The two-fermion dynamic quantities may exhibit van Hove singularities 
along the line uii'^k) = E^^^ki, k) where ki satisfies the equation 

J-i?(^)(fci,A;) = 0. (98) 
oki 

If for the solution of Eq. dMl) we also have d'^ E^3\ki,k) / dkl ^ the 
two-fermion dynamic structure factor shows the well-known square-root 
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Fig. 20. Two-fermion excitation continua (j = 1 (a, b), j = 2 (c, d), j = 3 (e, f)) for 
J = -1, 7 = 0.5, D = l,n = 0.5. Left panels: T ^ oo; right panels: T = 0. 

singularity. However, it may happen that d'^E^^^ki,k)/dkl — but 
d^E'^^\ki,k)/dk'l ^ 0. Then a van Hove singularity is characterized by 
the exponent 2/3. That is really the case, for example, for J = 1, 7 = 0.5, 
D = 1, f7 = 0.5 for j = 2 at fc = 1.0784. . .. We have dE^'^\ki,k) / dk^ = 
d^E^^^{ki,k)/dkl = 0, d^E^^\ki,k)/dkl ^ at fci = fcj = 2.1648.... 
Therefore in the e- vicinity of lo — 0.7859 . . . the two-fermion dynamic struc- 
ture factor should be proportional to |e|~^/'^. We mention that the singu- 
larity with this exponent is also present for D = C" - . 

In the zero-temperature case the effect of the Fermi functions involved 
in the C^-'^-functions (see Eq. ([96|l ) becomes important (right panels in 
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Fig. [^0]) . As a result the region of possible values of ki is contracted. Further 
details can be found in Ref. " . 

Finally, we mention the role of the iJ'^^^-functions (see Eq. (|96p ) for 
the two-fermion dynamic structure factors which are responsible for some 
specific features of the transverse dynamic structure factor (compare panels 
b, d, f in Fig. [20] and panel c in Fig. [T9|. 

We pass to the xx and yy dynamic structure factors (see panels a and 
b in Fig. I19p . These dynamic structure factors are many-fermion dynamic 
quantities and although they are not restricted to some region in the k-u) 
plane, they are rather small outside the two-fermion excitation continua 
(compare panels a and b with panel c in Fig. I19p . In the low-temperature 
limit the xx and yy dynamic structure factors show several washed-out exci- 
tation branches which are in correspondence with characteristic lines of the 
two-fermion excitation continua. In the high-temperature limit Sxx{k,uj) 
and Syy{k,oj) become fc-independent. 

It should be stressed that the constant frequency or constant wave- 
vector scans of the dynamic structure factors clearly manifest the presence 
of the Dzyaloshinskii-Moriya interaction and some easily recognized fea- 
tures of these quantities can be used for determining the Dzyaloshinskii- 
Moriya interaction. 

7. Square-Lattice Spin-1/2 Isotropic XY Model 

Let us discuss what kind of results for spin models can be obtained in two 
dimensions after applying the Jordan- Wigner transformation (Sec. [3]). We 
consider the spin-1/2 isotropic XY model on a spatially anisotropic square 
lattice with the Hamiltonian 



where J and J_l are the XX exchange interactions in the horizontal and 
vertical directions. Our aim is to calculate the transverse dynamic structure 
factor 




(99) 




(100) 
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We notice that the transverse dynamic structure factor SzzO^,^^) (|100p for 
the spin model ((99|) is related to the densit-^ensity dynamic structure 
factor of hard-core bosons on a square- lattice 

We apply the two-dimensional Jordan- Wigner transformation ((2T|) , ([23]) 
to the spin Hamiltonian (|99p . Moreover, we adopt the mean-field approach 
for the phase factors and change the gauge leaving the flux per elementary 
plaquette $o to be equal to tt. As a result we arrive at the Hamiltonian like 
(123, i.e. 



i=0 j=0 

+^ {dlAi+i - d^.A,+i)) ■ (101) 

The Hamiltonian (jlOll) contains the correct results in the one-dimensional 
limit when either J± = or J ~ (in the former case to recover the one- 
dimensional Hamiltonian (|16|1 (with = 0) one has to perform in addition 
a gauge transformation dj ^ = exp {mipi) fj ^ ipo = 0, — ^pi + i). 

The Hamiltonian poip can be diagonalized by performing 1) the Fourier 
transformation, dij = {l/ ^jN^Ny) J2k^,ky ^xp (i (f'xi + kyj)) dk^^ky, ka = 
27rn„/iV„, n„ = -iV„/2 + 1, . . T,No,/2 - 1, a = x,y, = Ny = 

^/N ^ CX3 is even, which yields 

iJ = i ^ |£^k| (cos 7k (blby: - a^aA + isin7k f^kflk - aj^^k 



l^'kl ^ y Jl cos^ fcj^ -I- J2 sin^ kx, 

J± COS ky . -^sinfe^r 
cos7k = ^r^, -n7k = ^^(102) 

with &k = dk^.ky and Ok = dk^±Tr,ky±Tr and 2) the Bogolyubov transforma- 
tion, ak = cos (7k/2) 6k + isin (7k/2)ak, /3k = sin (7k/2) 6k-i cos (7k/2) Ck, 
which yields 



Ak (aj^ak - /3i/3k) , 



Ak = l^^kl y ./i cos2 fcj^ + J2 sin2 fc:r > (103) 

(the prime near the sum in Eq. p03p means that k varies in the thermo- 
dynamic limit in the region — tt < ky < tt, —it + \ky\ < kx < tt — \ky\). 

The calculation of the transverse dynamic structure factor repeats the 
steps elaborated in some detail for the one-dimensional case. First, we use 
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the Wick-Bloch-dc Dominicis theorem to obtain the zz time-dependent spin 
correlation function 



2N ■ 

k 



((cos^ |i - i (_!)"+" _ 1) cos Y sin |i + sin^ |i 

•rik exp (iAkt) 

+ (sin^ |i + i (-1)"+™ - 1) cos |i sin |i + (-1)?'+? cos^ |i) 

• (1 - "k) cxp (-iAkt)) , 



IN ■ 

k 



• (fcos^ ^ + i (-1)"+™ - 1) COS ^ Sin ^ + sin2 

■ (1 - rik) exp (-iAkO 

+ (sin^ ^ - i (-1)"+'" ((-1)^+^ - 1) COS f sin |i + (-1)^+^ cos^ |) 

•nk exp (iAk0i)Q4) 

Then we plug Eq. (|104p into Eq. (jlOOp to obtain the following expression 
for the zz dynamic structure factor 

dfci„ dfcia; 



277 



cos2 ^'^^+'^^ ^'^^ Ki (1 - nki+k) (5 (c^ + Aki - Aki+k) 
+ (1 - J^ki) J^ki+ki^ (w - Aki + Aki+k)) 

I • 2 7ki+k ^ 7ki / /' I A I A ^ 

+ sm ^ (r^kinki+kC) (w + Aki + Aki+k) 

+ (1 - nkj (1 - riki+k) ^ii^- Aki - Aki+k))) ■ (105) 

One can easily convince oneself that Eq. (|105p contains the correct result in 
the one-dimensional limit (|35|) (with = 0) when either = or J = 0. 
In the two-dimensional case Eq. (jlOSp is an approximate formula for the 
transverse dynamic structure factor of the spin-1/2 isotropic XY model on 
a spatially anisotropic square lattice. 

Let us discuss the dynamic quantity obtained in some detail^. In 
Fig. [21] we show gray-scale plots for Szz^-,'-^) and in Fig. [22] we show fre- 
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Fig. 21. The zz dynamic structure factor Szz{kx, 0, lj) (gray-scale plots) for the square- 
lattice s = 1/2 XX model (|99j as it follows from Eq. I I105I I at T = (left column) and 
at T = 10 (right column). J = -1, Jx = -0.1 (a, b), Jx = -0.5 (c, d). 



2 



□ 




Fig. 22. Frequency dependence of the zz dynamic structure factor (I105I I for momentum 
transfer along the chain fc^; = 7r/2 as the interchain interaction changes (J = —1, = 
—0.1 (solid curves), J = —1, Jj^ = —0.5 (dashed curves), J = —1, J± = —0.9 (dotted 
curves)) at zero temperature T = (a) and high temperature T = 10 (b). 

quency profiles of SzzO^, ^) for a representative set of parameters. Formula 
(|105p implies the interpretation of iS'zz(k,w) as a two-fermion excitation 
quantity. As can be seen from Figs. [HI [221 5*2 2 (k, uj) exhibits several washed- 
out excitation branches which can be generated by two spinless fermions 
in accordance with (jlOSp . We begin with the low-temperature limit when 
only the fourth term in Eq. (jlOSp (which contains (1 — riki) (1 ~ 'T-ki-t-k)) 
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survives. Consider, for example, two fermions with ki = (0, 7r/2) — k/2 and 
ki + k = (0, 7r/2) + k/2 with the energy of the pair 



c^k = 2^ J2 sin^ y + Ji sin^ ^ (106) 

or two fermions with ki = (7r/2, 0) — k/2 and ki + k = (7i'/2, 0) + k/2 with 
the energy of the pair 



iU^ = 2-^/ J2 cos2 Y + Ji cos2 ^. (107) 



These modes are the well-known spin wavea^ clearly present at low tem- 
peratures (panels a and c in Fig. I2ip . Further, one can recognize the high 
frequency modes [ki = —k/2, ki = (7r/2, 7r/2) — k/2] with the dispersion 
relations 



c^k = 2-^/ J2 sin^ y + Ji cos2 ^, (108) 



co^ = 2^J^cos^^ + Jlsm^^. (109) 

Another set of high-frequency modes [ki = 0, ki = (7r/2,7r/2)] have the 
dispersion relations 



uJk^ Jj_ + \ ,P sin^ + J\ cos2 ky, (110) 



Wk = ./ + Y cos2 A;^ + ,/2 sin^ ky. (Ill) 
The low-frequency mode [ki = (0,7r/2)] with the dispersion relation 



= \J J'^ SVC? kx + J\svc? ky (112) 

(it is composed of two fermions, the energy of one of which equals zero) 
forms the low-frequency cutoff at zero temperature. Comparing left and 
right panels in Figs. [21] and [22] one can also see the modes which become 
visible only as temperature increases (at zero temperature they are for- 
bidden because of the Fermi factors in Eq. pOSp ). Putting kix = -^k^, 
kiy = tt/2 — kx for ky = and kix = ky, kiy = tt /2 — ky for k^ — we get 



uJo,k, = {\J^ + Ji-J]\smky\. (113) 
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This excitation branch contains most of the spectral weight at high tem- 
peratures (see panel d in Fig. [21] and panel b in Fig. [22]) . 

The established modes (|106p - (jll3|) manifest themselves as peaks, cusps 
or cutoffs in the frequency or wave- vector profiles of Szz (k, uj) . The fre- 
quency profiles depicted in Fig. [22] may be almost symmetric or asymmet- 
ric, they may resemble J-peaks or result from two coalesced peaks, they 
may gradually disappear or may be abruptly cut off. 

It is worthwhile to mention here some experimental studies on dynamic 
properties of two-dimensional quantum spin models, in particular, the neu- 
tron scattering experiments on Cs2CuCl4^^ (for a theor y o f dynamic cor- 
relations in the spin-liquid phase in CS2CUCI4 see Ref. l^^. CS2CUCI4 is 
a two-dimensional low-exchange quantum magnet. It has a layered crys- 
tal structure; in each layer the exchange paths form a triangle lattice with 
nonequivalent interactions along chains J = 0.374(5) meV and along zig-zag 
bonds J' = 0.34(3) J. The interlayer coupling is small J" = 0.045(5) J and it 
stabilizes the long-range magnetic order below Tn = 0.62(1) K. The neutron 
scattering measurements in the spin-liquid phase (i.e. above Tn but below 
J, J' when the two-dimensional magnetic layers are decoupled) clearly in- 
dicate that the dynamic correlations are dominated by highly dispersive 
excitation continua which is a characteristic signature of fractionalization 
of spin-1 spin waves into pairs of dcconfined spin- 1/2 spinous. Linear spin- 
wave theory including one- and two-magnonjarocesses cannot describe the 
continuum scattering. The proposed theories^ are based either on a quasi- 
one-dimensional approach (that immediately introduces spinon language) 
or on the explicitly two-dimensional resonating-valence-bond picture. 

As a final remark we recall that Eq. (I105|) contains the exact result (|35|) 
in the one-dimensional limit. On the other hand, Eq. (I105P gives an ap- 
proximate result in the two-dimensional case because of the mean-field de- 
scription of the phase factors which arise after fermionization. The adopted 
mean-field treatment neglects a complicated interaction between spinless 
fcrmions. In the case of the XXZ Heisenberg model the interaction be- 
tween fermions is present even within the adopted mean-field procedure 
due to the interaction between z spin components. The quartic terms in 
the fermionic Hamiltonian may be treated after making further approxima- 
tion (see references cited in Ref. and also Ref. . 
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8. Conclusions 



The Jordan- Wigner transformation which reahzes a spin-to-fermion map- 
ping was suggested as a rigorous framework for the description of spin-1/2 
XY chains in the early 1960s. In general, the Jordan- Wigner fermioniza- 
tion permits to map a system of interacting spins s — 1/2 onto a system 
of spinless fermions. It may happen that the spinless fermions are nonin- 
teracting. In this case this approach reveals an exactly solvable spin model. 
However, even for exactly solvable spin models not all 'simple' quantities of 
interest in spin language remain simple in fermionic language. For example, 
the z spin component attached to the site j, s|, becomes the product of 
two Fermi operators attached to this site, CjCj — 1/2. In contrast, the local 
spin operators sj, s|, become nonlocal objects in fermionic description 
involving a string of sites l,2,...,j (see Eqs. ([T2|) . ([TS])). This leads to 
some complications in studying the dynamics of fluctuations of these op- 
erators: the dynamics of fluctuations of operators which seem to be rather 
simple in spin language may be governed by many-particle correlations in 
fermionic language. As we have discussed in sections 21 [3 El the Jordan- 
Wigner fermionization approach permits to establish a number of rigorous 
results for the dynamics of spin-1/2 XY chains. Especially easy are the cases 
of two- and four-fcrmion dynamic quantities which are amenable mostly to 
analytical calculations. The case of many-fermion dynamic quantities is 
more complicated, however, these quantities can be examined numerically 
at very high precision. 

For more realistic spin-1/2 XX Z Heisenberg chains the Jordan- Wigner 
fermionization approach leads to a system of interacting spinless fermions. 
The simplest way to proceed in this case is to apply Hartree-Fock-like 

approximation^!^. If we are interested in low-energy physics only it might 

1 Q -| q on 

be helpful to apply the bosonization approach^" 

The results for one-dimensional quantum spin systems obtained using 
the Jordan- Wigner fermionization can be compared with the outcomes of al- 
ternative approaches: field-theoretic bosonization techniqueJI^IHIin! valid 
in the low-energy limit (see Fig. [72]) . Bethe ansatz calculations (for calcula- 
tion of dynamic structure factors of spin- 1 /2 XX Z chains see Refs. IMEOl) 
or exact diagonalization computations which, however, are restricted to 
small finite systems. 

For two-dimensional quantum spin models achievements are rather mod- 
est. In this case the Jordan- Wigner fermionization approach can provide 
an approximate theory; the simplest one treats in the mean-field spirit the 
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phase factors which arise after fermionization. 

We end up with a brief comment about the experimental relevance 
of some of the dynamic quantities calculated. They may be used for in- 
terpretation of the energy absorption in electron spin resonance (ESR) 
experimentJ^. Consider an ESR experiment in which the static magnetic 
field is directed along the z axis and the electromagnetic wave with the 
polarization in a _L z direction (say a — x) are applied to a magnetic sys- 
tem which is described as a spin-1/2 XY chain (ESR experiment in the 
standard Faraday configuration) . In such an ESR experiment one measures 
the intensity of the radiation absorption I{lu) as a function of frequency 
(jj > of the electromagnetic wave. Within the linear response theory the 
absorption intensity is written as 

/(w) CXC^9Xaa(0,Oj), (114) 

where '^Xaa{0, w) is the imaginary part of the aa component of the dynamic 
susceptibility Xaa{k,U!) at zero wave- vector fc = 0. We notice that 

3Xaa(0, io) - l^£^^|izM5^^(0, CO), (115) 

where the dynamic structure factor is defined by Eq. Thus, the pe- 
culiarities of the dynamic structure factor SaaiO,'^) caused, e.g., by the 
XY exchange interaction anisotropy, Dzyaloshinskii-Moriya interaction or 
dimerization should manifest themselves in ESR experiments. The time- 
dependent spin correlation functions taken at the same site or at the neigh- 
boring sites manifest themselves in the spin-lattice relaxation rate 1/Ti 
measured by nuclear magnetic resonance 

The activity in the field of the Jordan- Wigner fermionization approach 
has expanded much over the last few decades. Despite some limitations, 
the Jordan- Wigner fermionization approach has a wide range of applica- 
bility. Particularly attractive is that it allows one to handle complicated 
problems of low-dimensional quantum spin systems armed with relatively 
simple tools. It thus seems quite likely that it will continue to be used 
successfully in the coming years. 



Acknowledgments 

The author would like to thank T. Krokhmalskii, T. Verkholyak, J. Stolze, 
G. Miiller and H. Biittner in collaboration with whom the study of the 
dynamics was performed. He is grateful to T. Krokhmalskii for preparing 



27, 2008 16:39 



WSPC/Trim Size: 9in x 6in for Review Volume 



odcrzh'070612 



50 O. Derzhko 

all figures for the paper, many interesting conversations and helpful com- 
ments and suggestions. He thanks J. Stolze and T. Verkholyak for a criti- 
cal reading of the manuscript. NATO support is acknowledged (the grant 
reference number CBP.NUKR.CLG 982540, project "Dynamic Probes of 
Low-Dimensional Quantum Magnets" ) . The author acknowledges kind hos- 
pitality of the Organizers of the 43rd Karpacz Winter School of Theoretical 
Physics "Condensed Matter Physics in the Prime of XXI Century: Phe- 
nomena, Materials, Ideas, Methods" in L§,dek Zdroj in February 2007. 

References 

1. E. Ising, Z. Phys. 31, 253 (1925). 

2. L. Onsager, Phys. Rev. 65, 117 (1944). 

3. R. J. Baxter, Exactly Solved Models m Statistical Mechanics (Academic 
Press, London, 1982). 

4. H. Bethe, Z. Phys. 71, 205 (1931). 

5. M. Karbach and G. Miiller, arXivicond-ma t/9809162[ 



M. Karbach, K. Hu and G. Miiller, arXiv:cond-mat/9809163! 
M. Karbach, K. Hu and G. Miiller, arXiv:cond-mat/0008018 

6. M. Kenzelmann, R. Coldea, D. A. Tennant, D. Visser, M. Hofmann, 
P. Smeibidl and Z. Tylczynski, Phys. Rev. B65, 144432 (2002). 

7. J.-S. Caux, F. H. L. Essler and U. Low, Phys. Rev. B68, 134431 (2003); 
D. V. Dmitriev and V. Ya. Krivnov, Phys. Rev. B70, 144414 (2004); 

D. V. Dmitriev and V. Ya. Krivnov, Pis'ma Zh. Eksp. Tear. Fiz. 80, 349 
(2004). 

8. M. D'lorio, R. L. Armstrong and D. R. Taylor, Phys. Rev. B27, 1664 (1983); 
M. D'lorio, U. Glaus and E. StoU, Solid State Commun. 47, 313 (1983). 

9. P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

10. Th. J. Siskens and P. Mazur, Physica A71, 560 (1974). 

11. E. Lieb, T. Schultz and D. Mattis, Ann. Phys. (N.Y.) 16, 407 (1961). 

12. S. Katsura, Phys. Rev. 127, 1508 (1962); 
S. Katsura, Phys. Rev. 129, 2835 (1963). 

13. V. M. Kontorovich and V. M. Tsukernik, Sov. Phys. JETP 25, 960 (1967); 
Th. J. Siskens, H. W. Capel and K. J. F. Gaemers, Physica A79, 259 (1975). 

14. M. Suzuki, Prog. Theor. Phys. 46, 1337 (1971); 

D. Gottlieb and J. Rossler, Phys. Rev. B60, 9232 (1999); 

0. Derzhko, J. Richter and V. Derzhko, Ann. Phys. (Leipzig) 8, SL49 (1999); 

1. Titvinidze and G. I. Japaridze, Eur. Phys. J. B32, 383 (2003). 

15. J. P. de Lima and L. L. Gongalves, J. Magn. Magn. Mater. 206, 135 (1999); 
O. Derzhko, J. Richter and O. Zaburannyi, Physica A285, 495 (2000); 

O. Derzhko, J. Richter, T. Krokhmalskii and O. Zaburannyi, Phys. Rev. 
E69, 066112 (2004); 

J. P. de Lima, L. L. Gongalves and T. F. A. Alves, Phys. Rev. B75, 214406 
(2007). 

16. H. Nishimori, Phys. Lett. AlOO, 239 (1984); 



27, 2008 16:39 



WSPC/Trim Size: 9in x 6in for Review Volume 



odcrzh'070612 



Jordan- Wigner fermionization 51 



O. Derzhko and J. Richter, Phys. Rev. B55, 14298 (1997); 
O. Derzhko and J. Richter, Phys. Rev. B59, 100 (1999). 

17. L. N. Bulaevskii, Zh. Eksp. Theor. Fiz. 44, 1008 (1963) [Sov. Phys. JETP 
17, 684 (1963)]. 

18. A. Luther and I. Peschel, Phys. Rev. B12, 3908 (1975). 

19. H. J. Schulz, Phys. Rev. B34, 6372 (1986). 

20. I. Affleck, Field theory methods and quantum critical phenomena. In: Fields, 
Strings and Critical Phenomena, Ed. E. Brezin and J. Zinn- Justin (Elsevier, 
Amsterdam, 1989) pp. 563-640; 

J. von Delft and H. Schoeller, Ann. Phys. (Leipzig) 7, 225 (1998); 
S. Rao and D. Sen, arXiv:cond-mat/0005492 

21. H. H. Fu, K. L. Yao and Z. L. Liu, Phys. Rev. B73, 104454 (2006). 

22. E. Fradkin, Phys. Rev. Lett. 63, 322 (1989); 

E. Fradkin, Field Theories of Condensed Matter Systems (Addison-Wesley, 
1991). 

23. Y. R. Wang, Phys. Rev. B43, 3786 (1991); 
Y. R. Wang, Phys. Rev. B43, 13774 (1991); 
Y. R. Wang, Phys. Rev. B46, 151 (1992). 

24. M. Azzouz, Phys. Rev. B48, 6136 (1993); 

M. Azzouz, L. Chen and S. Moukouri, Phys. Rev. B50, 6233 (1994); 

M. Azzouz and C. Bourbonnais, Phys. Rev. B53, 5090 (1996); 

M. Azzouz, B. Dumoulin and A. Benyoussef, Phys. Rev. B55, R11957 

(1997). 

25. D. C. Cabra and G. L. Rossini, Phys. Rev. B69, 184425 (2004); 

C. A. Lames, D. C. Cabra, M. D. Grynberg and G. L. Rossini, Phys. Rev. 
B74, 224435 (2006). 

26. C. D.Batista and G. Ortiz, Phys. Rev. Lett. 86, 1082 (2001). 

27. S. V. Dobrov, J. Phys. A36, L503 (2003). 

28. M. N. Kiselev, D. N. Aristov and K. Kikoin, Phys. Rev. B71, 092404 (2005). 

29. O. Derzhko, Journal of Physical Studies (L'viv) 5, 49 (2001). 

30. G. Misguich, Th. Jolicoeur and S. M. Girvin, Phys. Rev. Lett. 87, 097203 
(2001); 

D. C. Cabra, M. D. Grynberg, R C. W. Holdsworth and P. Pujol, Phys. Rev. 
B65, 094418 (2002); 

M.-C. Chang and M.-F. Yang, Phys. Rev. B66, 184416 (2002). 

31. A. Kitaev, Ann. Phys. (N.Y.) 321, 2 (2006). 

32. X.-Y. Feng, G.-M. Zhang and T. Xiang, Phys. Rev. Lett. 98, 087204 (2007). 

33. Th. Niemejer, Physica 36, 377 (1967). 

34. G. Miiller, H. Thomas, H. Beck and J. Bonner, Phys. Rev. B24, 1429 (1981); 

G. Miiller, H. Thomas, M. W. Puga and H. Beck, J. Phys. C14, 3399 (1981). 

35. J. H. Taylor and G. Miiller, Physica A130, 1 (1985). 

36. J. Lorenzana and G. A. Sawatzky, Phys. Rev. Lett. 74, 1867 (1995); 
J. Lorenzana and G. A. Sawatzky, Phys. Rev. B52, 9576 (1995); 

H. Suzuura, H. Yasuhara, A. Furusaki, N. Nagaosa and Y. Tokura, 
Phys. Rev. Lett. 76, 2579 (1996); 

Yongmin Yu, G. Miiller and V. S. Viswanath, Phys. Rev. B54, 9242 (1996); 



27, 2008 16:39 



WSPC/Trim Size: 9in x 6in for Review Volume 



oderzh'070612 



52 O. Derzhko 



J. Lorenzana and R. Eder, Phys. Rev. B55, R3358 (1997); 
R. Werner, Phys. Rev. B63, 174416 (2001). 

37. H. G. Vaidya and C. A. Tracy, Phystca A92, 1 (1978). 

38. O. Derzhko, T. Krokhrnalskii, J. Stolze and G. Miillcr, Phys. Rev. B71, 
104432 (2005); 

O. Derzhko, T. Krokhrnalskii, J. Stolze and G. Miiller, Physica B378-380, 
445 (2006). 

39. T. Barnes, Phys. Rev. B67, 024412 (2003). 

40. U. Brandt and K. Jacoby, Z. Phys. B25, 181 (1976); 
U. Brandt and K. Jacoby, Z. Phys. B26, 245 (1977). 

41. H. W. Capel and J. H. H. Perk, Physica A87, 211 (1977); 
J. H. H. Perk and H. W. Capel, Physica AlOO, 1 (1980). 

42. H. B. Cruz and L. L. Gongalves, J. Phys. CIA, 2785 (1981). 

43. B. M. McCoy, E. Barouch and D. B. Abraham, Phys. Rev. A4, 2331 (1971); 
J. Lajzerowicz and P. Pfeuty, Phys. Rev. Bll, 4560 (1975); 

B. M. McCoy, J. H. H. Perk and R. E. Shrock, Nucl. Phys. B220, 35 (1983); 
B. M. McCoy, J. H. H. Perk and R. E. Shrock, Nucl. Phys. B220, 269 (1983); 
G. Miiller and R. E. Shrock, Phys. Rev. B29, 288 (1984). 

44. F. Colomo, A. G. Izergin, V. E. Korepin and V. Tognetti, Phys. Lett. A169, 
243 (1992); 

A. R. Its, A. G. Izergin, V. E. Korepin and N. A. Slavnov, Phys. Rev. Lett. 

70, 1704 (1993); 

A. R. Its, A. G. Izergin, V. E. Korepin and N. A. Slavnov, Phys. Rev. Lett. 
70, 2357 (1993); 

F. Colomo, A. G. Izergin and V. Tognetti, J. Phys. A30, 361 (1997). 

45. H. Asakawa, Physica A233, 39 (1996). 

46. O. Derzhko and T. Krokhrnalskii, Phys. Rev. B56, 11659 (1997); 

O. Derzhko and T. Krokhrnalskii, phys. stat. sol. (h) 208, 221 (1998). 

47. A. P. Young and H. Rieger, Phys. Rev. B53, 8486 (1996); 

S. Sachdev and A. P. Young, Phys. Rev. Lett. 78, 2220 (1997); 
A. P. Young, Phys. Rev. B56, 11691 (1997). 

48. J. Stolze, A. Noppert and G. Miiller, Phys. Rev. B52, 4319 (1995). 

49. J. Stolze and M. Vogel, Phys. Rev. B61, 4026 (2000). 

50. Y. Maeda and M. Oshikawa, Phys. Rev. B67, 224424 (2003). 

51. X. Jia and S. Chakravarty, Phys. Rev. B74, 172414 (2006). 

52. O. Derzhko, T. Krokhmalskii and J. Stolze, J. Phys. A33, 3063 (2000). 

53. K. Fabricius, U. Low and J. Stolze, Phys. Rev. B55, 5833 (1997). 

54. O. Derzhko, T. Krokhmalskii and J. Stolze, J. Phys. A35, 3573 (2002). 

55. K. Kawasaki, N. Maya, A. Kouzuki and K. Nakamura, J. Phys. Soc. Jpn. 
66, 839 (1997). 

56. A. Kouzuki, K. Kawasaki and K. Nakamura, Phys. Rev. B60, 12874 (1999). 

57. W. Yu and S. Haas, Phys. Rev. B62, 344 (2000). 

58. M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 (1997); 
D. N. Aristov and S. V. Maleyev, Phys. Rev. B62, R751 (2000); 

O. Derzhko, J. Richter and O. Zaburannyi, J. Phys.: Condens. Matter 12, 
8661 (2000). 



27, 2008 16:39 



WSPC/Trim Size: 9in x 6in for Review Volume 



odcrzh'070612 



Jordan- Wigner fermionization 53 

59. O. Derzhko and T. Verkholyak J. Phys. Soc. Jpn. 75, 104711 (2006). 

60. M. Oshikawa and I. Affleck, Phys. Rev. B65, 134410 (2002). 

61. J. Sirker, Phys. Rev. B73, 224424 (2006). 

62. O. Derzhko, T. Verkholyak, T. Krokhmalskii and H. Biittner, Phys. Rev. 
B73, 214407 (2006); 

O. Derzhko, T. Verkholyak, T. Krokhmalskii and H. Biittner, Physica 
B378-380, 443 (2006). 

63. T. Matsubara and H. Matsuda, Prog. Theor. Phys. 16, 569 (1956); 
H. Matsuda and T. Matsubara, Prog. Theor. Phys. 17, 19 (1957); 
H. Matsuda, Prog. Theor. Phys. 18, 357 (1957); 

T. Morita, Prog. Theor. Phys. 18, 462 (1957). 

64. O. Derzhko and T. Krokhmalskii, Physica B337, 357 (2003); 

O. Derzhko and T. Krokhmalskii, Acta Physica Polonica B32, 3421 (2001); 
O. Derzhko and T. Krokhmalskii, J. Magn. Magn. Mater. 232-245, 778 
(2002); 

O. Derzhko and T. Krokhmalskii, Czechoslovak Journal of Physics 55, 601 
(2005). 

65. G. Gomez-Santos and J. D. Joannopoulos, Phys. Rev. B36, 8707 (1987). 

66. R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, B. Dorner and 
Z. Tylczynski, Phys. Rev. Lett. 79, 151 (1997); 

R. Coldea, D. A. Tennant, A. M. Tsvelik and Z. Tylczynski, Phys. Rev. Lett. 
86, 1335 (2001); 

R. Coldea, D. A. Tennant, K. Habicht, P. Smeibidl, C. Wolters and Z. Tyl- 
czynski, Phys. Rev. Lett. 88, 137203 (2002); 

R. Coldea, D. A. Tennant and Z. Tylczynski, Phys. Rev. B68, 134424 (2003). 

67. M. Bocquet, F. H. L. Essler, A. M. Tsvehk and A. O. Gogolin, Phys. Rev. 
B64, 094425 (2001); 

M. Bocquet, Phys. Rev. B65, 184415 (2002); 

C. H. Chung, J. B. Marston and R. H. McKenzie, J. Phys.: Condens. Matter 
13, 5159 (2001); 

C.-H. Chung, K. Voelker and Y. B. K im, Phys. Rev. B68, 094412 (2003); 
Yi Zhou and X.-G. Wen, 'arXiv:cond-mat/0210662 

68. A. Lopez, A. G. Rojo and E. Fradkin, Phys. Rev. B49, 15139 (1994). 

69. M. Karbach, G. Miiller, A. H. Bougourzi, A. Fledderjohann and K.- 
H. Mutter, Phys. Rev. B55, 12510 (1997); 

M. Karbach, D. Biegel and G. Miiller, Phys. Rev. B66, 054405 (2002). 

70. J.-S. Caux and J. M. Maillet, Phys. Rev. Lett. 95, 077201 (2005); 

J.-S. Caux, R. Hagemans and J. M. Maillet, J. Stat. Mech.: Theor. Exp., 
P09003 (2005); 

J.-S. Caux and R. Hagemans, J. Stat. Mech.: Theor. Exp., P12013 (2006). 



